Ultra-cold Fermi gases with resonant dipole-dipole interaction 
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The superfluid phases in the resonant dipolar Fermi gases are investigated by the standard mean- 
field theory In contrast to the crossover from Bose-Einstein condensation (BEC) to Bardeen-Cooper- 
Schrieffer (BCS) superfluid in the Fermi gases with the isotropic interactions, the resonant dipolar 
interaction leads to two completely different BEC phases of the tight-binding Fermi molecules on 
both sides of the resonance, which are characterized by two order parameters with the distinct 
internal symmetries. We point that near the resonance, the two competitive phases can coexist, and 
an emergent relative phase between the two order parameters spontaneously breaks the time-reversal 
symmetry, which could be observed in the momentum resolved rf-spectroscopy. 
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Introduction. — The unprecedented experimental pro- 
gresses p]-[5| in creating quantum gases of fermions with 
large dipole moment has stimulated extensive investiga- 
tions of dipolar Fermi gases. Owing to the long-range and 
anisotropic nature of the dipole-dipole interaction (DDI), 
new quantum phenomena emerge in the dipolar Fermi 
gases, e.g., the ferro-nematic phase [6], the novel static 
and dynamical properties in the normal phase [7], and 
the p-wave dominated BCS superfluids induced by the 
partially attractive DDI (sl-Hlj. Of particular interest, 
the recent studies demonstrated that in the two-species 
dipolar Fermi gases the competition between the short- 
range contact interaction and DDI led to the coexistence 
of singlet- and triplet-paired superfluids [l2l-[l4j. 

The above theoretical work focus on the Fermi gases 
with DDI in the weak coupling regime, for which the 
scattering amplitude is replaced by the first Born ap- 
proximation of the bare DDI potential [l5j|. Recent ex- 
perimental progresses toward creating ultracold gases of 
polar molecules |l|-l3| and Rydberg atoms [4] provide an 
opportunity for exploring the dipolar effects in the strong 
coupling regime. Previous studies on the low-energy scat- 
tering between two polarized dipoles indicate that, due 
to the shape resonances, the scattering lengths diverge 
at certain strengths of the DDI |sL [l6l4l8| . In particular, 
these shape resonances may take place simultaneously in 
multiple scattering channels as the DDI couples partial 
waves with different angular momenta [I?). This multi- 
channel resonance (MCR), described by a matrix of the 
scattering lengths, is in striking contrast to the Feshbach 
resonance (see, e.g., Ref. [19| and references therein). 

In this Letter, we investigate the novel superfluid 
phases in a two-species dipolar Fermi gas with resonant 
DDI. Based on the analysis to the low-energy scatter- 
ing of two polarized dipoles, we propose an effective 
two-body interaction potential and a model to describe 
the MCR. It is found that a two-body bound state is 
formed on either side of the resonance. Consequently, 



for a many-body system, the Bose-Einstein condensates 
of the tight-binding molecules dominates on both sides of 
the resonance. Since these two molecular states possess 
distinct internal symmetries, the system experiences a 
phase transition across the MCR. Moreover, near the res- 
onance, the two competitive phases coexist and a relative 
phase between the two order parameters emerges, which 
spontaneously breaks the time-reversal symmetry (TRS). 
Similar mixed order parameters were studied in the high- 
T c superconductors |2Q|-[24j. Finally, we study the quasi- 
particle spectral function of the system in order to ex- 
plore the possibility of the experimental detection using 
the momentum-resolved (MR) rf-spectroscopy [25|, [26| . 

Model. — We consider an ultracold gas of two-species 
dipolar fermions in a box of volume V. For simplicity, 
we assume that = = n with n a being the number 
density of the spin-a particle. The total Hamiltonian 
of the system can be decomposed into H = Ho + Hi, 
where, in the momentum space, Ho = a e^c^c^a 
with Cka being the annihilation operator of the spin-a 
fermion and = k 2 / (2M) — \i with M being the mass of 
the particle and \i the chemical potential. Furthermore, 
the interaction Hamiltonian takes the form 

a/3 kpq 

where U(k) = 47r<i 2 (3cos 2 #k — l)/3 is the Fourier trans- 
form of the bare DDI between two polarized (along z- 
axis) dipoles with d being the dipole moment and #k the 
angle between k and the z-axis. We note that the short- 
range interaction can be straightforwardly included in the 
two-body interaction potential U [13|, [l5| . 

Effective interaction potential. — In order to obtain 
valid results in the strong coupling regime, the bare DDI 
has to be renormalized. To this end, we consider the 
scattering between two polarized dipoles. In terms of the 
K- matrix, /C(fc), the scattering amplitude can be gener- 
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ally expressed as 



f(k',k)\ k=k , = 4tt ( r=TL 
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where the subscripts and superscripts of the bracket de- 
note the matrix elements. Specifically, for DDI, since 
U oc Y20, K>\™ is nonzero only when \l — l'\ < 2. In ad- 
dition, because the DDI conserves the projection of the 
angular momentum, /C is diagonal with respect to the 
magnetic quantum number m. Of particular importance, 
in the zero energy limit, the matrix elements JC\™(k)/k 
are all nonvanishing [8|, LL6| , which results in finite scat- 
tering lengths affi 
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To proceed further, we construct a matrix A whose 
elements are defined by the scattering lengths as Affi = 
i l '~ l a\™\ Assuming that the matrix A is diagonalized in 
the orthonormal basis Wj m (k) = djiYimfe) with the 
corresponding eigenvalues Xj m . In this new basis, the 
scattering amplitude Eq. (|2]) can be reexpressed as 
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where fj m = — l/(A - r ^ + ik). It is clear that the eigen- 
value Xjm represents the effective scattering length in the 
scattering channel Wj m (k.). We can now construct a sep- 
arable effective interaction potential as 



U'(k,k>) = 4nJ2g jm w jm (k)w* m (k'), 
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where the coupling constants g jm satisfy the renormal- 
ization condition 
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The correctness of the effective potential U' can be ver- 
ified as it reproduces the two-body scattering amplitude 
Eq. (|2]), therefore, it is valid in both weak and strong 
coupling regimes. We point out that two other pseudo- 
potentials applicable to the strong DDI were also pro- 
posed previously (27], [28| . 

Two-body physics in the resonance regime. - With 
the effective potential U' ', it can be shown that if A jm 
is positive, the scattering channel Wj m (k) may support 
a two-body bound state with binding energy E^jm = 
— l/(MA| m ) (see, e.g. Ref [29]). Indeed, numerical cal- 
culations show that as one tunes the dipole moment d of 
the colliding particles, the scattering lengths affi appear 
to have various resonances BKEa. Those shape res- 
onances indicate the DDI indeed supports bound states. 
Noting that, within the same / and V manifold, affi with 
m = is always the largest term [8], here and henceforth, 
we restrict our analysis to m = and drop the index m. 



To be specific, we consider the collision between a spin- 
t and a spin-^ particles in the spin-singlet channel. In 
general, the shape resonances induced by the DDI are 
well-separated [17j , which allows us to focus on a partic- 
ular resonance, say the resonance of aoo at dipole mo- 
ment d r . Interestingly, Kanjilal and Blume found that 
a resonance on ao2 also takes place at the same position 
d r [17j. In principle, because the DDI couples the partial 
waves \l — 1'\ < 2, resonances should occur on all aw with 
l,l f — even. Therefore, unlike the Feshbach resonance, 
the shape resonance induced by the DDI may occur si- 
multaneously in multiple scattering channels. Utilizing 
the fact that the widths of the resonances decrease with 
increasing I + V [17|, we propose a minimal model for the 
MCR by assuming that it is described by the matrix 
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where all other scattering lengths are assumed to be zero. 
We note that the exact behavior of ^4 s d depends on the 
details of the short-range physics of the colliding parti- 
cles. However, as will be shown, some general properties 
for the MCR can be obtained by analyzing this simplest 
model. 

The matrix *4 s d can be easily diago nalized to y ield the 
eigenvalues Ai, 2 = [&oo ± sgn(a 2) \/ a oo + 4a 02]/ 2 anc ^ 
the corresponding eigenstates (k) = [si 5 2^oo(k) + 



^20(k)]/ 
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1, where si ?2 = — (y zb \/y 2 + 4)/2 with 



y = «oo/«02- Clearly, independent of the values of the 
scattering lengths, only one of the eigenvalues can be 
positive, i.e., there only exists a single bound state for 
any given set of (aoo,«o2)- More specifically, on the up- 
per (<%) , C&02 ) P arame ^ er plane, we have Ai > 0, which 
leads to the binding energy Eb = —l/MXf and the 
angular distribution of the bound state wave function 
I i^i (k)| 2 . On the other hand, A2 is positive on the lower 
(<2oo 5 a o2 1 ) plane. Consequently, the angular distribution 
of the bound state wave function becomes \w2 (k) | 2 • Since 
w 2 (k) is orthogonal to wi(k), the bound state changes its 
microscopic symmetry when ao2 changes its sign. One 
can also carry out a similar analysis when aoo changes 
its sign. Therefore, we expect that a quantum phase 
transition arises as the system crosses an MCR. 

We remark that the above discussion can be easily gen- 
eralized to the MCR between two spin polarized par- 
ticles by replacing the even orbital angular momentum 
quantum numbers with odd ones. In this case, the sim- 
plest model for the MCR is characterized by the matrix 



A 



-ai3 



ajl q 13 I. The properties of the two-body 



bound state can be obtained straightforwardly. 

Many-body physics across an MCR. — Now, we turn to 
investigate the many-particle properties of an ultra-cold 
gas of two-species fermionic dipoles within the regime 
of an MCR. At zero temperature, the gap and phase 
fluctuations can be ignored. The system can then be 
treated by using the standard BCS mean-field theory. 
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FIG. 1. (Color Online) Order parameters (upper row) and 
chemical potential (lower row) for a two-species dipolar Fermi 
gas in the resonance regime with y = — 1 (left column) and 1 
(right column). Here, sf — k F /(2M) is the Fermi energy. 



To start, we define the order parameters for the chan- 
nels Wj(k) as Aj = Airgj J2k w j(M) ( c -k|Ck^) /V, which 
satisfy the coupled gap equations 
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|A(k)| 2 is the quasi-particle excita- 
tion energy with A(k) = J2j ^j w jO^) being the order 
parameter. Furthermore, to completely determine the 
gaps Aj, one needs the density equation, which, for the 
spin-balanced system, takes the form 



3^ 



d 3 p 
(2tt)3 



(7) 



where kF = (67r 2 n) 1//3 is the Fermi momentum. Equa- 
tions (j6j) and (0 form a closed set of equations for the 
order parameters Aj and the chemical potential fi. Even 
though, the order parameters Aj are generally complex 
numbers, it can be shown that only their absolute values 
|Aj| and the relative phase <j) = arg(Ai) — arg(A2) are 
relevant. 

For a given set of scattering lengths aoo and ao2, the 
self-consistent Eqs. (0 and ([7]) can be solved numerically. 
In general, there exist multiple solutions corresponding 
to different local minima of the energy density, 
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The true ground state can be identified by comparing 6 
corresponding to the different solutions. 

In principle, to obtain a full picture of the superfluidity 
in the resonant dipolar gases, one should numerically find 
the gaps for the scattering lengths covering the entire 
(clqo,clq2) plane. Here, for simplicity, we fix the ratio 
V — &00/&02 of the scattering lengths and treat a^ 1 as a 



free parameter. The assumption that y is fixed implies 
that the widths of the resonances in aoo and ao2 are the 
same, which seems rather unrealistic, as the resonance in 
ao2 is generally more narrow than that in aoo- However, 
the physics present here stands as long as the widths of 
the two resonances become comparable. 

In Fig. [TJ we present the order parameters and the 
chemical potential of the system across an MCR for 
y = ±1. As can be seen, away from the resonance, the 
condensation of the molecular state [Ai (A 2 ) for ao2 > 
(< 0)] always dominates. From the angular distributions 
of the order parameter A(k), the Cooper pairs shows 
completely different internal symmetries on the different 
sides of the resonance. In this sense, a phase transition 
takes place as the system crosses the resonance. This re- 
sult is in contrast to the Feshbach resonance, for which 
only a smooth crossover is experienced |3(J Ell. Near the 
resonance, a^ 1 ~ 0, the amplitudes of Ai and A2 become 
comparable such that the two order parameters coexist 
in the system. Of particular interesting, a nonzero rela- 
tive phase also develops in this region. As a result, A(k) 
becomes complex, which indicates that the TRS is spon- 
taneously broken. In fact, such mixed order parameters 
with nonzero relative phases have been extensively stud- 
ied in the high-T c superconductors [20| . 

Further evidence about the phase transition is provided 
by the chemical potential which becomes negative away 
from the resonance. This indicates that strongly cou- 
pled BECs form on the both sides of the resonance. In 
addition, \i reaches its maximal value when the system 
approaches the resonanance, which explicitly shows the 
competition between these two phases. As a comparison, 
we point out that for the BCS-BEC crossover induced by 
the Feshbach resonance, the chemical potential monoton- 
ically decreases from the Fermi energy to a negative value 
as the system passes through the resonant regime from 
the BCS side to the BEC side. 

It is also worthwhile to note that the sign of aoo de- 
termines the detailed behavior of the order parameter for 
the BCS state. For negative (positive) aoo, the amplitude 
of the BCS state quickly (smoothly) decays to near zero 
when the system moves away from the resonance. 

Replacing the matrix *4 s d by *4 p f, one can generalize 
the above analysis to a single component dipolar Fermi 
gas in the MCR regime. Indeed, it is found that other 
than the angular distribution of the order parameter, the 
behaviors of the gaps and the chemical potential are very 
similar to those of a two-component gas. 

Experimental detection. — Now, we explore the possi- 
bility for the experimental detection of the fermionic su- 
perfluid with MCR. A potential technique for detecting 
the phase with mixed order parameters and the associ- 
ated phase transition is the MR rf-spectroscopy [25|. In 
such experiment, an rf pulse drives the transition from 
one of the two spin states to an unoccupied third spin 
state, and the photoemission spectroscopy measures the 
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FIG. 2. (Color Online) Quasi-particle spectral function, 
lnA(k, uj), at cos#k = ±0.39 (upper row) and ±0.9 (lower 
row) for a superfluid in the multiple-channel resonant regime. 
From the left to the right columns, l/(fcpao2) = — 1, 0, and 
1, respectively. Other parameters are y = 1 and j/sf = 0.2. 








I cos k | 



FIG. 3. (Color Online) y dependence of the angular distribu- 
tions for |wi(0k)| (left panel) and \w2{0\^)\ (right panel). 



quasi-particle spectral function 

u 2 

A(k,u>) = - V t ^ (8) 

V J 7T (cJ-s£ k ) 2 + 7 2 ' V J 

where ± = (1 ± £^/E\^)/2 and 7 is the energy resolu- 
tion. 

As an example, we present lnA(k, oS) in Fig. [2] for the 
ratio of the scattering lengths y = 1. At the left end 
of the resonance {l/hp^m ~ —1), the system is domi- 
nated by the order parameter A2^2(k). It can be shown 
that the gap is closed at cos#£ ~ ±0.39 and reaches its 
maximum value at cos#£* ~ ±0.9. This properties is 
clearly demonstrated in the left column of Fig. [2j On 



the other hand, for l/kpa^ ~ 1 5 the order parameter 
roughly becomes Aiiui(k), for which the gap is closed 
at and opens up at The behavior of the spectral 
function can then be used to identify the phase transi- 
tion. Moreover, as shown in the middle column of Fig. [2j 
the gap presents for both and in the resonance 
regime (l/kpCLm ~ 0). In fact, in this regime, it can be 
shown that the gap is nonzero for arbitrary 0^. 

For an arbitrary y, it is still possible to distinguish 
the two phases away from the resonance by using the 
spectral function. To demonstrate this, we plot the y 
dependence of the angular distributions of the order pa- 
rameters Wj(k) in Fig. [3] It can be shown that for 
V < Vi - 1.79 (y > y 2 - -0.22), w x (w 2 ) always 
has a zero at |cos#k| > 0.58 (< 0.58). Therefore, for 
2/2 < V < the spectral function possesses similar fea- 
tures as those displayed in Fig. [2j While for y < y 2 and 
y > Vi j only at one end of the resonance one can observe 
that the gap is closed. The angle #k at which the gap 
vanishes can be used to identify the phase. 

Finally, we note that the similar analysis for the spec- 
tral function of the MCR with odd angular momenta can 
be carried out straightforwardly, for which the gap always 
vanishes at cos#k = 0. 

Conclusions. — We have studied the superfluid phases 
of a two-species dipolar fermionic gas across an MCR. 
Based on the analysis to the low energy scattering of two 
polarized dipoles, we derive a separable effective potential 
for the DDI which is valid in the strong coupling regime. 
Subsequently, we proposed a minimal model which de- 
scribes the MCR. We then studied the two-body and the 
many-body physics across an MCR. It was found that 
condensates of the molecular states with distinct internal 
symmetries form on both sides of the resonance, indicat- 
ing that a phase transition arises when the system passes 
through the resonance. Near the resonance, the two com- 
petitive order parameters coexist. In addition, a relative 
phase between those two orders emerges, which spon- 
taneously breaks the TRS. We also studied the quasi- 
particle spectral function of the system in oder to ex- 
plore the possibility of the experimental detection using 
the MR rf-spectroscopy. Finally, we point out that due 
to the spontaneous TRS breaking, the resonant dipolar 
Fermi gas provide an opportunity for studying the Ma- 
jorana edge modes in the spin-polarized dipolar Fermi 
gases [32[. 
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